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I. INTRODUCTION 

One of the unsolved problems of the Standard Model(SM) is the hierarchy problem. The little Higgs models 1, 
HI H, 0] are introduced to solve the hierarchy problem by stabilizing the Higgs mass by a collective symmetry 
breaking mechanism due to the cancelation of divergent loops by appearance of new particles as a consequence of 
extra symmetries. The phenomenology of the little Higgs models are widely discussed in literature (for reviews 
see [1, [1, 0])j an d constraints on little Higgs model parameters are studied d, [T3, HH, E3, EH- The little Higgs 
models are also expected to give new significant signatures in future high energy colliders and studied in references 
. Il7|. due to the new particles which are predicted by these models. Also the Zl associated production 
of SM Higgs boson at e + e~ colliders in the framework of the littlest Higgs model is studied in [TH, 

In this work we examined the production of single and doubly ch arg ed scalars associated with Zl boson at future 
linear e + e~ colliders, namely, International Linear Collider (ILC) [2(| and Compact Linear Collider (CLIC) [2l| in 
the context of littlest Higgs modelQ. We examined the dependance of total and differential cross sections of the 
processes to the littlest Higgs model parameters at the range allowed by electroweak presicion measurements. Also 
we analyzed the decays of scalars in the context of lepton flavor violation. We found that the production rates of 
the single charged scalar pairs associated with Zl are less than the Zl associated production of doubly charged 
, scalars, but both channels will be achieved at e + e~ colliders at yfs > 2TeV. 

The layout of this paper is as follows: In section 2 first we give a brief review of the littlest Higgs model and then 
we calculate the pair production of charged scalars at e + e~ colliders. In this section we also calculate the total and 
differential cross sections of the productions of charged scalars with Zl- Section 3 contains our numerical results 
and discussions. 



II. THEORETICAL FRAMEWORK 



Before examining the pair production of charged scalars with Zl, we introduce a few words on the littlest Higgs 
model based on reference [6]. In the littlest Higgs model global symmetry SU(5) is broken spontaneously to SO(5) 
at an energy scale / ~ ITeV leaving 14 nambu goldstone bosons(NGB) corresponding to broken symmetries 
represented by the goldstone boson matrix of non linear sigma model (nlsm) such as: 
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Il(x) = J2u a (x)X a (1) 

a=l 
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where X a are broken generators of SU(5). The vacuum bases triggering the symmetry breaking can be chosen as 




and the E field is defined as 

E(s) = e m /fZ a e m /f = e 2ffl /'£ , 
and the effective lagrangian of the E field is: 

f 2 

£ E = ^-Tr|P u E| 2 . 



(2) 
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In the littlest Higgs model SU(5) contains the gauged subgroup [SU(2) 1 <g> E/(l)i] <g> [SU(2) 2 <g> U(l) 2 ] of SU(5), 
defining the covariant derivative such as: 



(5) 



»=i 



where and W^i are the gauge fields, g[ and <?, are the corresponding couplings of U(l)i and SU(2)i respectively. 
The vev of E field in the lagrangian breaks (SU(2) (8 U(l)) 2 symmetry to diagonal subgroup (SU{2) ® U(l)) of 
SM. As a consequence symmetry breaking, gauge bosons gain mass by eating the four of the NGBs. The mixing 
angles between the SU(2) subgroups and between the U(l) subgroups are defined respectively as: 



(-12 



92 



(6) 



The usual electroweak symmetry breaking occurs by the vev of the Higgs potential written by Coleman Weinberg 
method for scalars. By EWSB vector bosons get extra mixings due to vacuum expectation values of h doublet 

2 

and 4> triplet. Again by diagonalizing the mass matrices, the final masses of vector bosons to the order of jj are 
expressed as@|: 
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where m w = gv/2, m z 
weak mixing angles: 



/ // 2 '2_i_ 2 /2\ 

gv/(2c w ) and = %gg' ^ & i C ,i ,1 ± i\ - In these equations s w and c w are the usual 
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The parameters u and u' in equation are the vacuum expectation values of scalar doublet and triplet given as; 
(h°) = v/y2 where v — 2A6GeV and {icjP) — v' < |j, bounded by electroweak precision data. Also diagonalizing 

the mass matrix for scalars the physical states are found to be the SM Higgs scalar H, the neutral scalar the 
neutral pseudo scalar P , and the charged scalars <j) + and cf> ++ . The masses of these scalars are degenerate, and in 
terms of Higgs mass can be expressed as: 
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The scalar fermion interactions in the model are written in yukawa lagrangian preserving gauge symmetries of the 
model for SM leptons and quarks, including the third generation having an extra singlet, the T quark. In this 
work, leptons are charged under both U(l) groups, with corresponding hypercharges of Y\ and Y 2 . The restriction 
for Y\ and Y% is that Y\ + Y2 should reproduce U (l)y hypercharge Y of SM, thus Y\ = xY and Y% = (1 — xjY can 
be written. Due to gauge invariance, x can be taken as 3 / 5[ a Il3||. Also for light fermions, a Majorano type mass 
term can be implemented in yukawa lagrangian H3, HH HEf which results in lepton flavor violation by unit two, 
such as: 

C LFV = iY ZJ Lf 0C-%-+h.c, (10) 

where Li are the lepton doublets ( I Vi), and Yij is the yukawa coupling with Yu = Y and Y^^y) = Y' . The 
values of yukawa couplings Y and Y' are restricted by the current constraints on the neutrino masses (26|. given 
as; Mij = YijV 1 ~ 10~ 10 GeV[22|]. Since the vacuum expectation value v' has only an upper bound; v' < lGeV, Yij 
can be taken up to order of unity without making v' unnaturally small. 

In the littlest Higgs model the symmetry breaking scale /, and the mixing angles s, s' are not restricted by the 
model. These parameters are constrained by observables of electroweak presicion data[l, 0, [HI HH, EH, For 
the values of the symmetry breaking scale ITeV < f < 2TeV, mixing angles are in the range 0.75 < s < 0.99 and 
0.6 < s' < 0.75, for 2TeV < f < 3TeV they have acceptable values in the range 0.6 < s < 0.99 and 0.6 < s' < 0.8, 
for 3TeV < f < 4TeV they are in the range 0.4 < s < 0.99 and 0.6 < s' < 0.85, and for the higher values 
of the symmetry breaking scale, i.e. / > ATeV , the mixing angles are less restricted and they are in the range 
0.15 < s < 0.99 and 0.4 < s' < 0.9 0. 

After these preliminary remarks lets start to calculate the cross sections of the pair productions for the charged 
scalars with at e + e~ colliders. In the littlest Higgs model there are four neutral vector bosons; the SM Z boson 
Zl, massless photon Al, and two new bosons Zh and Ah- Their couplings to fermions are written as i^i^igvi + 
9Ailb) where i = 1, 2, 3, 4 corresponds to Zl, Zh, Ah and Al respectively. The couplings of gauge vector to electron 
positron pairs are given in table HI where e = ^/Aira, y e = | for anomaly cancelations, x^ = — sc(c 2 — s 2 ) 

and x§ = —ji-s'c'(c' 2 — s 12 ). It is seen that vector and axial vector couplings of SM Z^e + e~ vertex also gets 
contributions from littlest Higgs model. As a result total decay widths of SM vector bosons also gets corrections of 
the order ji, since the decay widths of vectors to fermion couples are written as; T(Vi — > //) = ^(ffy + g^^Vi 
where TV = 3 for quarks, and N = 1 for fermions. The total decay widths of the new vectors are given as (l7| : 
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The new scalars and pseudo scalars also contribute to the analysis done in this study. Since these new scalars 

have lepton flavor violating decay modes, their total widths will depend on the yukawa couplings Yu = Y and 

Yij^ijtj) = Y' if the flavor violating term in the yukawa lagrangian is considered. The decay widths of scalars are 
given as[22l|: 

?/ 2 A/f 3 ^ 

' V| 2 M + ^|r'| 2 M , (12) 
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where we neglected the decays into quarks since they are very small compared to other terms. The first term 
in the decay width proportional to \Y\ 2 corresponds to decays into same family of leptons and the second term 
proportional to \Y'\ 2 correspond to the decays into different families. The branching ratios of scalars decaying to 
same family of leptons is denoted as BR[Y] and to leptons of different flavor is BR[Y']. In this work the values of 
the yukawa mixings are taken to be 10~ 4 < Y < 1 and Y' < 10~ 4 . 

The Feynman diagrams for the Zl associated scalar pair production processes are given in figure [1] The first 
ten of the Feynman diagrams in figure [1] are for the both processes e + e~ — > ZL<j) + <f>~ and e + e~ — > Z^0 ++ </) . 
The last four ones only contribute to the associated production of single charged scalars. Now we present the 
amplitudes for these Feynman diagrams. The relevant Feynman rules for the vertices for single charged scalars 
are given in tables [Til Mil and IIVI and for the doubly charged scalars in tables [V] and IVII The amplitudes for the 
associated production of doubly charged scalars are get by replacing the couplings with the primed ones. In the 
calculations of cross sections, terms are expanded up to order of v 2 / f 2 . 

The first four Feynman diagrams contributing to both processes contain only the e + e~Vi vertices. The amplitudes 
for the first two diagrams in figure [TJ are written as: 



Mi = u[-p 2 ]ij IJ ,(g Vl + g Al j 5 )e^[p 3 ] t -Yi^g Vi u[p 1 }(i)^--iEf' p '(p 4 -ps)o, (13) 
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TABLE I: The vector and axial vector couplings of ee with vector bosons. Feynman rules for eeVi vertices are given as 



i/j 


vertices 




1/1 


4> + 4>~ ZlZl 




2/1 


4> + (f>~ ZhZl 


0(^ 2 // 2 ) ~ 


3/1 


<j> + (t>-AHZ L 


o„-gg' (o rJ -s'") „2 
ZJ c w 2s'c' 


4/1 


<f> + 4>~ AlZl 





TABLE II: Feynman rules for the four point 4> + 4> ViVj vertices 0. 
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W "' iEf<>>(p 4 -p 5 ) a , (14) 



g' 2 - Af 2 + iM^T 

where g = p% — P3, (f = Pa +P5 and coefficients Ef are given in table If f f \ E'^ for doubly charged scalars given in 
table rvT|) . The amplitudes for diagrams 3 and 4 in figure UJ are written as: 



M 3 = u[-p 2 ]ij fJ ,(gvJ-^-iEt' t> (p 4 - p 5 )^ a (gv 1 + 9A 1 ls)e a [p?,}u[p 1 ], (15) 
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iAij a (gvi + 9Ai75h a \p3\u\pi\, (16) 

<r 

where q = p% — pa, q' = p 4 + p$ and coefficients E^{E'f ) are given in table UTTljVip . The rest of the amplitudes 
from 5 to 14 in figure UJ are s channel processes with q = p±+p2- For the diagrams 5 to 8 in figure [TJ corresponding 
amplitudes 5 to 8 have a heavy scalar as a propagator, and proportional to the square of the couplings of two 
vectors with a charged scalar. These amplitudes are written as: 

M 5 = u[-p 2 ]i'y fJ ,g Vi u\p 1 ]i—iE%' t> (p 4 _ q ') v _ iEf^i-q' - p 3 ,) a e a [p 3 }, (17) 
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TABLE III: Feynman rules for 4> + ij> Vi vertices^. 
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TABLE IV: The Feynman rules for (f> + ViW '• vertices. Their couplings arc given in the form iBf- g^v where g^ v carries the 
lorentz indices of vectors and j = 1, 2 denotes Wl, Wh respectively [g|. 



i/j 


vertices 




1/1 


cj> ++ 4>—z L z L 


2i4(l-24) 2 5M „ 


2/1 


<f> ++ <f,~-Z H ZL 




3/1 


<j> ++ <j } —A H Z L 


-2^ l£ ^P(l-24)^ 


4/1 


<t> ++ <i>—A L Z L 


4*e£(l - 24) 5mi , 



TABLE V: Feynman rules for the four point 



'ViVj vertices 0| . 



gW i 
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Ms = ^2u[~p2}ijf,{g Vi + 9Ail5)u[pi}'i 
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where Ef^(E'^) are given in table IIIlt|VI|) . and q' = q — p 5 for amplitudes 6 and 8 and q' = q — Pi for 5 and 7. The 
diagrams 9 and 10 in figure [1] are sub processes in which the four point couplings of vectors and scalars contribute. 
Their amplitudes are given as: 



Mg = u[-p 2 ]ij fl g Vi u[pi]i-- I -iC£{'g ua e' T \p3], 



(21) 
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i/j 


vertices 
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TABLE VI: Feynman rules for 4> ++ 4> V t vertices [i 




FIG. 1: Feynman diagrams contributing to e + e — > Zl4> + 4> an d e+e ~~ v Zl4> ++( I ) m littlest Higgs model. The last 
four diagrams only contribute to the single charged scalars. 

where Cf-(C' i j) are four point couplings given in table ITTlfVf . The amplitudes 11 to 14 corresponds to diagrams 
11 to 14 in figure[T]are only for single charged scalars, where an electric charge is carried by a W propagator. These 
sub processes are proportional to the couplings of a vector with a charged scalar and a charged vector Wj , where 
j = 1,2 corresponds to Wl and Wh respectively. Their amplitudes are given as: 
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where Bf- are given in table ITVl and q' = q — p§ for amplitudes 11 and 13 and q' = q — pi for 12 and 14. 

The numerical calculations of cross sections of the pair production of charged scalars are performed by CalcHep 
generator [28T ]. 

III. RESULTS AND DISCUSSIONS 

In this section we present and discuss our results for the Zl associated pair production of charged scalars. In 
performing the numerical calculations, we take the electromagnetic coupling constant e = V Ana = 0.092, the Higgs 
mass Mh = l20GeV and the mass of the standard model bosons Mz L = 9lGeV, Mw L — 80GeV and the SM 
mixing angle sw = 0.47 using the recent data[27j|. In the calculations, we ignored v 2 / f 2 terms in the couplings, 
since we are not dealing with the corrections to a SM process. 

For the pair production of single charged heavy scalars associated with Zl , the differential cross sections 
versus energy of the Zl boson graphs for fixed values of mixing angle parameters s/s' at symmetry breaking scale 
/ = ITeV at total center of mass energy y/s = 3TeV appropriate for CLIC are presented in figure [2]^ a) . The total 
cross sections for these processes for different fixed values of the parameters s/s' — 0.8/0.6, 0.8/0.7, 0.95/0.6, 0.5/0.1 
are presented in table I VIII The dependance of the total cross section on yfs for these processes are presented in 
figure[3l The differential cross section gets its maximum value about 0.1^^? when s/s' — 0.5/0.1, correspondingly 
for total cross section we obtain that 59/6. This means that thousands of productions per year at high integrated 
luminosity of 100/6" 1 are expected. But this parameter set is not allowed when symmetry breaking scale / = ITeV 
by electroweak precision data. For parameters s/s' = 0.8/0.6, 0.8/0.7, 0.95/0.6, the peak values of differential cross 



e + , e" -> Z L ,* + ,<J>" e + , e" -> Z L ,4> ++ ,*" 
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FIG. 2: Differential cross section vs. Ez graphs for processes e + e~ — + Zz,4> (j>~ (a) and e + e~ — > Zl4> ++ 4> (b) for the fixed 
value of parameters s/s' : 0.8/0.6, 0.8/0.7, 0.95/0.6, 0.5/0.1 at / = lOOOGeV at s/s = 3TeV. 
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0.5/0.1 
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TABLE VII: The total cross sections in fb for pair production of charged scalars associated with Zl for / = ITeV and at 
= 3TeV. 

sections are obtained at the order of 10~ 4 ^7 for low Ez values, Ez ~ 100GeV. The total cross section is calculated 
as 0.04/6. This result implies 1 ~ 10 events per year accessible for a collider luminosity of 100/6 -1 . The single 
charged scalar <p~ has leptonic decay modes Uvi and Uvj with branching ratios proportional to \Y 2 \ and \Y' 2 \ 
respectively. It also has decay modes to SM particles; W^H, W^Zl pairs. The SM decays dominates when the 
yukawa coupling Y is either zero or very small (Y <C 1). For higher values of yukawa coupling (Y ~ 1) the leptonic 
modes (Uvi) of </> + _ couple dominates the final states, giving a signature of missing energy of the neutrino anti 
neutrino couple, and a lepton anti lepton pair of same or different flavor plus Zl- These signals can be lepton 
flavor violating but the observation is challenging due to the low production rates of the single charged pair. 

For the associated pair production of doubly charged scalars within Zl the differential cross sections with respect 
to Ez arc plotted in figure ED^b) for different fixed values of mixing angles at \/s — 3TeV, at symmetry breaking 
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FIG. 3: Total cross section vs. \fs graphs for Zl associated pair production of charged scalars at / = ITeV when parameters 
s/s' : 0.8/0.7. 
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FIG. 4: Total cross section vs. ^fs graphs for Zl associated pair production of charged scalars at (a)/ = 2.5TeV , and at 
(b)/ = 3.5TeV when parameters s/s' : 0.8/0.7. 



scale / = ITeV. The dependence of total cross section on the center of mass energy of this production process 
is plotted in figure O and the numerical values of total cross sections are presented in table I VIII for parameters 
of interest. The differential cross section of the production process reaches its maximum when model parameters 
s/s' = 0.5/0.1 about O.l^^r for Ez ~ lOOGey, resulting a total cross section of 84/6. This will give about 
8000 events per year for high luminosities such as 100/6 -1 . But this remarkable event numbers are out of reach, 
since s/s' = 0.5/0.1 is out favored by electroweak presicion data for / = ITeV. For the electroweak allowed 
parameters s/s' = 0.8/0.6,0.8/0.7,0.95/0.6 at / = ITeV at y/s — 3TeV, the differential cross section gets low 
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FIG. 5: Total cross section vs. y/s graphs for Zl associated pair production of charged scalars at / = 5TeV when parameters 
(a)s/s' : 0.8/0.7 and (b)s/s' : 0.6/0.4. 
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values at the order of lO -4 ^^. The resulting cross sections are calculated by integrating over Ez, and found about 
0.4 ~ 0.8 x fb (table EH) resulting 40 ~ 80 events per year for integrated luminosity of lOOfb^ 1 . For yfs < 2TeV, 
this production channel is not reachable. In the littlest Higgs model, <j) ++ has decays to charged vectors W^W^ 
and also to leptons iflj proportional to squares of the values of the yukawa couplings; \Y 2 \ for same families and 
|F' 2 | for different lepton families when lepton violating modes are considered. So this channel provides final signals 
for doubly charged scalar discovery and lepton flavor violation. The final states of the doubly charged scalar pairs 
dominantly contain leptonic modes hhlflj, semi leptonic modes lihW^W^ and to standard model charged vector 

pair W^W^W£W£ depending on Y and Y', while Zl dominantly decays to jets carrying the energy at the order 
of masses of the scalar pair. For / = ITeV the leptonic branching ratio of doubly charged scalars can reach values 
close to 1 for Y — ► 1, independent from Y'. If the value of the yukawa coupling Y is high enough (Y ~ 1), the 
number of final state lepton flavor violating signals such as; Z^lililjlj , can reach up to 50 events per year for 
luminosities of 100/6 -1 , which can be directly detectable free from backgrounds. 

Finally, we have also analyzed the behavior of the production processes for higher values of /. The dependence 
total cross section on y/s when s/s' — 0.8/0.7 for / = 2.5TeV,3.5TeV,5TeV are plotted in figures Ufa), BJb) and 
13 a), respectively. It is seen that the maximum value of total cross sections, production rates and final lepton 
flavor violating signals for these / values for both processes remain same, but the required energy is shifted to 
y/s = 3.8TeV(5.5TeV){7.5TeV} for / = 2.5TeV(3.5TeV){5TeV}, due to the increase in heavy scalar masses. For 
/ = 5TeV, we have also analyzed the case s/s' = 0.6/0.4, since the parameters are less constrained. The dependence 
of total cross sections for associated production of both single and doubly charged pairs are plotted in figure \5j[h) ■ 
It is seen that the total cross sections are increased by order one. For the single charged pair, the production cross 
section is calculated as 0.9/6 for yfs > lOTeV, resulting 90 production events for luminosities 100/6 -1 . For the 
double charged scalar pair, the total cross section is 5 ~ 9fb for yfs > lOTeV, giving 500 ~ 900 productions at 
luminosities 100/6 -1 . In this case the final number of four lepton signals (lililflj) will be around 600 for higher 
values of yukawa coupling (Y ~ 1). If the values of yukawa coupling is smaller, in the region 0.1 < Y < 0.3, the 
double charged pair will decay into semi leptonic modes, such as W^W^kh, resulting 200 ~ 400 signals accessible 
for luminosities of 100/6 -1 . 

In conclusion, we found that at an e + e~ collider of y/s > 2TeV with a luminosity of 100/6 -1 , the Zl associated 
production of charged scalars will be in the reach, being the single charged pair is quite challenging due to low 
production rates, and the production of doubly charged scalar pair more promising for the electroweak allowed 
parameters at / = ITeV. The final states will contain lepton flavor violating signals if the value of yukawa 
coupling Y is close to unity. For larger values of / the mixing angles s/s' are less constrained, e.g. for / = 5TeV 
and s/s' = 0.6/0.4, the production rates increases allowing remarkable final lepton number violating events for 
0.1 < Y < 1, but for these set of parameters the center of mass energy of the colliders should be increased. 
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